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ABSTRACT 

PSR B1951 +32 is a y-ray pulsar detected by the Energetic Gamma Ray Experiment Telescope (EGRET) 
and identified with the 39.5 ms radio pulsar in the supernova remnant CTB 80. The EGRET data 
shows no evidence for a spectral turnover. Here we report on the first observations of PSR B1951+32 
beyond 30 GeV. The observations were carried out with the 10m y-ray telescope at the Whipple Ob¬ 
servatory on Mt. Hopkins, Arizona. In 8.1 hours of observation we find no evidence for steady or 
periodic emission from PSR B1951+32 above ~ 260 GeV. PEux upper limits are derived and com¬ 
pared with model extrapolations from lower energies and the predictions of emission models. 

INTRODUCTION 

The pursuit of Very High Energy (VHE) astrophysics has resulted in the discovery of five sources, 
of which three are associated with young spin-powered pulsars. VHE emission has been detected 
from the direction of the Crab Nebula (Vacanti et.al., 1991), the Vela pulsar (Takanori 1996) and PSR 
B1706-44 (Kifune et ah, 1995) but no evidence has been found for periodic emission at these energies 
in these experiments. 

PSR B1951+32 has been detected as a pulsating X-ray source (Safi-Harb et ah, 1995) and as 
a high energy y-ray pulsar at E > 100 MeV at the radio period (Ramanamurthy et ah, 1995). It can 
be inferred from the five pulsars seen in the MeV to GeV y-ray region that longer period or older 
( ~ 10^ years ) pulsars have a greater fraction of spin down energy emitted as high energy y-rays. 
The best fit outer gap model of Zhang and Cheng (1997) suggests that PSR B1951+32 should emit 
detectable levels of TeV y-rays (Pigure|2l). The multiwavelength spectrum of PSR B1951+32 (Eigure 
lb) indicates a maximum power per decade at energies consistent with a few GeV and still rising at 
10 GeV. These factors make PSR B1951+32 a good candidate for observations with the ACT above 
100 GeV. 

OBSERVATIONS 

The observations of PSR B1951+32 reported here were acquired with the 10m reflector located at the 
Whipple Observatory on Mt. Hopkins in Arizona. A total of 14 Tracking runs and 4 On/Off pairs 
taken between 13th May, 1996 and 17th July, 1996 constitute the database for all subsequent discus¬ 
sion. The total On source observing time is 8.1 hrs. The radio position (J2000) of PSR B1951+32 (a 
= 19^ 52"^ 58.25'^, 5= 32° 52' 40.9") was assumed for the subsequent timing analysis. 


PSR 

P 

msec 

o 

5i 

VI 

V 

1 

Distance 

kpc 

logioB 

Gauss 

EogioE 

ergs/s 

B1951+32 

39.53 

5.8494 

2.5 

11.69 

36.57 


ANALYSIS AND RESULTS 
Standard Analysis 

The event selection criteria are collectively called Supercuts95 and a detailed description can be found 
elsewhere (Catanese et ah, 1995). Supercuts95 raises the effective energy threshold of the detector 
with its software trigger and size cuts. PSR B1951+32 appears to have a steep spectrum at EGRET 
energies and since the pulsar spectrum is expected to cut off, it behooves us to reduce the threshold 
of our analysis to search for a lower energy signal. The dominant background at lower energies is 
due to muons whose images appear in the camera as arcs and can be discriminated by a cut on their 
large length/size values. Hence the selection criteria used Supercuts95 on images with sizes larger 
than 400 p.e. and Smallcuts (Table|2l) for images with sizes less than 400 p.e. No steady emission is 

Table 2: Parameter ranges for selecting y-ray images 


Parameter 

Supercuts95 

Smallcuts 

length 

0°16-0°30 

unchanged 

width 

0°073-0°15 

unchanged 

distance 

0°51 - ri 

unchanged 

alpha 

< 15° 

unchanged 

maxi 

> 100 p.e. 

45 p.e. - 100 p.e. 

max2 

> 80 p.e. 

45 p.e. - 80 p.e. 

size 

> 400 p.e. 

0 - 400 p.e. 

length/size 

not used 

< 7.5 X 10-4 °/p.e. 


Table 3: Selected Events for Steady Emission analysis 


Selection 

Source Events 
a< 15° 

Background Events 
a< 15° 

Excess 

Significance 

Supercuts95 

292 

254 

38 

L16o 

Smallcuts 

618 

672 

-54 

-l.lOa 

Supercuts95 + Smallcuts 

910 

926 

-16 

-0.24O 


detected from PSR B1951+32 and 3a flux upper limits are displayed in TablelH The effective area for 
Supercuts95, that was used to calculate the upper limit, was taken as Agjp ~ 3.5 x 10* cm^: the same 
area was used for the dataset that resulted from a combination of Supercuts95 and Smallcuts although 
here there is more systematic uncertainty. The energy threshold was obtained from simulations and 
extrapolating the Crab Nebula y-ray rate for each set of cuts used assuming a spectrum ~ 

Periodic Analysis 

The arrival times of the Cerenkov events were registered by a GPS clock with an absolute resolution 
of 250 psQC. An oscillator calibrated by GPS second marks was used to interpolate to a resolution of 



















system barycenter and folded to produce the phases, of the events modulo the pulse period. The 
ephemeris frequency parameters used were v= 25.2963719901267 s^^ and v=-3.73940x lO^^^s^^, 
at the epoch to=JD 2450177.5. This frequency was extrapolated 72 days to obtain a timing solution 
relevant to the epoch of observation. The datasets, however, were taken within the validity interval of 
the above ephemeris. 

To check the Whipple Observatory timing systems an optical observation of the Crab pulsar 
was undertaken on December 2nd (UT) 1996 using the 10m reflector. The phase analysis of the event 
arrival times yielded a clear detection of optical Crab pulsar signal which is in phase with the radio 
pulse and demonstrates the validity of the timing, data acquisition and software in the presence of a 
pulsed signal. No evidence of pulsed emission from PSR B1951+32 at the radio period exists. To 
calculate a pulsed flux upper limit we assumed the same pulse profile as seen at EGRET energies, i.e. 
with the phase range for the main pulse and secondary pulse as 0.12 - 0.22 and 0.48 - 0.74 respectively 
(Ramanamurthy et ah, 1995). 

Table 4: Integral Elux Upper limits 



Steady Emission 

Periodic Emission 

Threshold 


cm^^s”^) 

(cm-^s-i) 

(GeV) 

Supercuts95 

0.97 X 10-" 

3.7x 10-1^ 

> 370 

Supercuts95 + Smallcuts 

1.95x 10-11 

6.7 X 10-12 

> 260 


DISCUSSION 

PSR B1951+32 is surrounded by a compact nebula 
which gives a plerionic nature to the supernova rem¬ 
nant, CTB80. X-ray plerions are good candidates for 
VHE emission since the electrons responsible for neb¬ 
ular synchrotron X-rays should also create VHE y-rays 
via the inverse Compton (IC) process. It is expected 
that for plerions, such as that associated with PSR 
B1951+32 where the density of nebular synchrotron 
photons is too low for SSC to take place, detectable 
VHE emission should be produced by the IC scatter¬ 
ing of the 2.7K cosmic microwave background by the 
same electrons radiating synchrotron X-ray photons. 
Interpreting the unpulsed X-ray emission form CTB80 
as the synchrotron emission from a plerion, the esti¬ 
mated IC flux > I TeV is 6.6 x 10^^^ TeV/cm^/s/TeV 
(De lager et. ah, 1995). This represents the lower limit 
on the IC flux since there can be other sources of soft 
photons in addition to the microwave background. 

To model the pulsed high energy spectrum, a 
function of the form 

dN^/dE = KE 



Log V (Hz) 


Eig. 1: The pulsed energy spectrum of PSR 
B1951+32. The Whipple limit is indicated 
as a filled square at 370 GeV. (See text for 
details). 


was used where E is the photon energy, E is the photon spectral index and Eq is the cut off energy. 
The pulsed upper limit reported here is two orders of magnitude lower than the extrapolated EGRET 
power law. Equation (2) was used to extrapolate the EGRET spectrum to VHE energies constrained 





emission (Figure la). 

The strength of the cut off provides a good dis¬ 
criminant between the various pulsar emission mod¬ 
els. The status of current observations and the derived 
cutoff discussed above indicates that the cutoff is be¬ 
yond 10 GeV. In polar cap models this would indicate 
a sharp cutoff since the pair production optical depth 
increases exponentially with photon energy (Harding 
1997). However, it is not possible to constrain the 
shape of the cut off with the non detection of pulsed 
TeV flux reported here. The most relavant compari¬ 
son of the Whipple upper limit with emission models 
is the outer gap model of Zhang and Cheng (see Fig¬ 
ure El). This model includes the effect of geometry in 
the treatment of pulsed emission via a parameter a = 
xiri, the radial distance to the synchrotron emitting re¬ 
gion near the outer gap, r, as a function of the light 
cylinder radius Our pulsed upper limits are consis¬ 
tent with the outer gap model if a > 0.6 implying an 
emission region far out in the magnetosphere. 

The result reported here is the first observation 
of PSR B1951-1-32 beyond 30 GeV. PSR B1951-1-32 
exhibits very similar spectral behavior and morpho¬ 
logical features, such as an associated synchrotron nebula, to PSR B1706-44 (Finley et ah, 1997). If 
these factors are any indication of similar emission mechanisms in pulsars then the lack of unpulsed 
emission from PSR B1951-1-32 is puzzling considering that PSR B1706-44 was detected as a VHE 
source of unpulsed emission > 1 TeV (Kifune et ah, 1995). Lack of pulsed emission indicates that the 
processes producing pulsed high energy photons over two decades of energy in the EGRET energy 
range somehow become ineffective over a decade of energy to result in a lack of VHE y-rays. The 
low magnetic field of PSR B1951-1-32 relative to the average pulsar field implies that attenuation of 
y-rays by magnetic absorption is not a likely explanation for the non-detection. 
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Fig. 2: Predicted pulsed y-ray flux of PSR 
B1951+32 from the Zhang and Cheng outer- 
gap model. The solid, dot-dash and dashed 
curves correspond to a=0.5, 0.6, 0.7 respec¬ 
tively. (See text for details). 



